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Finding control fields (pulse sequences) that can compensate for the dispersion in the parameters 
governing the evolution of a quantum system is an important problem in coherent spectroscopy 
and quantum information processing. The use of composite pulses for compensating dispersion in 
system dynamics is widely known and applied. In this paper, we make explicit the key aspects of 
the dynamics that makes such a compensation possible. We highlight the role of Lie algebras and 
non-commutativity in the design of a compensating pulse sequence. Finally we investigate three 
common dispersions in NMR spectroscopy, the Larmor dispersion, rf-inhomogeneity and strength 
of couplings between the spins. 
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I. INTRODUCTION 



Many applications in control of quantum systems in- 
volve controlling a large ensemble by using the same con- 
trol field. In practice, the elements of the ensemble could 
show variation in the parameters that govern the dynam- 
ics of the system. For example, in magnetic resonance 
experiments, the spins of an ensemble may have large 
dispersion in their natural frequencies (Larmor disper- 
sion), strength of applied rf- field (rf-inhomogeneity) and 
the relaxation rates of the spins. In solid state NMR 
spectroscopy of powders, the random distribution of ori- 
entations of inter-nuclear vectors of coupled spins within 
an ensemble leads to a distribution of coupling strengths 
|l0j . A canonical problem in control of quantum ensem- 
bles is to develop external excitations that can simulta- 
neously steer the ensemble of systems with variation in 
their internal parameters from an initial state to a de- 
sired final state U H H ■ These are called compen- 
sating pulse sequences as they can compensate for the 
dispersion in the system dynamics. From the standpoint 
of mathematical control theory, the challenge is to simul- 
taneously steer a continuum of systems between points 
of interest with the same control signal. Typical applica- 
tions are the design of excitation and inversion pulses in 
NMR spectroscopy in the presence of larmor dispersion 
and rf-inhomogeneity Q1I1I10jIEI1IH1] or the trans- 
fer of coherence or polarization in coupled spin ensemble 
with variations in the coupling strengths :9|. In many 
cases of practical interest, one wants to find a control 
field that prepares the final state as some desired func- 
tion of the parameter. For example, slice selective excita- 



tion and inversion pulses in magnetic resonance imaging 
[H El El The problem of desi gning excitations 
that can compensate for dispersion in the dynamics is a 
well studied subject in NMR spectroscopy and extensive 
literature exists on the subject of compositepulses that 
correct for dispersion in system dynamics Q,Hy,Q- The 
focus of this paper is not to construct a new compensat- 
ing pulse sequence but rather to highlight the aspects of 
system dynamics that make such a compensation possi- 
ble and give proofs of existence of a compensating pulse 
sequence. Our final goal is to understand what kind of 
dispersions can and cannot be corrected. 

To fix ideas, consider an ensemble of noninteracting 
spin i in a static field Bq along z axis and a transverse 
rf-field, {A(t) cos(</>(i)), A(t) sin(</>(i))), in the x — y plane. 
Let x, y, z represent the coordinates of the unit vector in 
direction of the net magnetization vector of the ensemble. 
The dispersion in the amplitude of the rf-field is given by 
a dispersion parameter e such that A(t) = eAo(t) where 
e G [1 — S, 1 + S] , for 5 > 0. Similarly there is dispersion in 
the larmor frequency to around a nominal value luq, i.e., 
lu — loo = Auj G [— B, B]. In a rotating frame rotating 
with frequency loq, the Bloch equations take the form 
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Consider now the problem of designing controls u(t) 
and v(t) that simultaneously steer an ensemble of such 
systems with dispersion in their natural frequency and 
strength of rf-field from an initial state (x, y, z) = (0, 0, 1) 
to a final state (x,y,z) = (1,0,0) |£g. This prob- 
lem raises interesting questions about controllability, i.e., 
showing that inspite of bounds on the strength of rf- 
field, ^/u 2 (t) + v 2 (t) < A max , there exist excitations 
(u(t),v(t)), which simultaneously steer all the systems 
with dispersion in Alu and e, to a ball of desired radius r 
around the final state (1, 0, 0) in a finite time (which may 
depend on A max , B, S, and r). These are control prob- 
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lems involving infinite dimensional systems with special 
structure. Besides steering the ensemble between two 
points, we can ask for a control that steers an initial 
distribution of the ensemble to a final distribution, i.e., 
if X(t) denote the units vector (x(t),y(t), z(t)), consider 
the problem of steering an initial distribution X(Au>, e, 0) 
to a target function X(Aw, e, T) by appropriate choice of 
controls in equation (Q. If a system with dispersion in 
parameters can be steered between states that have de- 
pendency on the dispersion parameter, then we say that 
the system is ensemble controllable with respect to these 
parameters. A more formal definition will appear later 
in the paper. 

This paper is organized as follows. In the following sec- 
tion, we introduce the key ideas and through examples, 
highlight the role of Lie brackets and non-commutativity 
in the design of a compensating control. In section 3, we 
show that the Bloch equations with bounded con- 
trols, u(t) and v(t) are ensemble controllable in the pres- 
ence of Larmor dispersion and rf-inhomogeneity. Finally 
in section 4, we investigate in some generality, the notion 
of ensemble controllability for linear control systems and 
a class of nonlinear control systems. 

II. LIE BRACKETS AND ENSEMBLE 
CONTROLLABILITY 

Example 1: Main Concept To fix ideas, we begin by 
considering Bloch equations with only rf-inhomogeneity 
and no Larmor dispersion. 

X = e(u(t)ft y + v(t)Q x )X 



where 
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are the generators of rotation around x, y and z axis, 

respectively. 

Observe for small dt, the evolution U\{\/~dt) — 

exp(— eily^/dt) exp(— efl x Vdt) exp(eCl y V~dt) enp(efl x Vdt) 

to leading order in dt is given by I+(dt)[ett y , efl x ], i.e., we 
can synthesize the generator [efl x ,efl y ] = e 2 fl z , by back 
and forth maneuver in the directly accessible directions 
Q x and Q y . 

Similarly, the leading order term in the evolution 
U 2 = Ui(-Vdl)eKp(-eflydt)U 1 (Vdt)exp(eQ y dt). 

is [efl y , [e£l x , e£ly]] = e 3 fl x . Therefore by successive Lie 
brackets, we can synthesize terms of the type e 2k+1 Q x . 
Now using {efl x ,e 3 il x , . . . ,e 2n+1 il x } as generators, we 
can produce an evolution 

n 

exp^cfe^aj, 

fc=0 



where n and the coefficient Cfc can be chosen so that 

n 
fc=0 

for all e £ [1 — S, 1 + 8] . Hence we can generate an evo- 
lution exp(8£l x ) for all e to any desired accuracy. There- 
fore, we achieve robustness with dispersion to e by gen- 
erating suitable Lie brackets. Similar arguments show 
that we can generate any evolution exp(/3£lj,) and as a 
result any three dimensional rotation in a robust way. It 
is also now easy to see that we can synthesize rotation 
with a desired functional dependency on the parameter e. 
Parametrize a rotation in £ SO (3) by the Eulcr angles 
{a, (3, 7) such that = exp(afi x ) exp(/3f2 g ) exp(7fi x ). 
Given continuous function (a(e), /3(e), 7(e)), of e, we can 
find polynomials that approximate a(e), /3(e) and 7(e) 
arbitrarily well and use these to generate a desired rota- 
tion 0(e) as a function of e. Hence there exists a control 
field that maps a smooth initial distribution X e (0) to a 
target distribution Xf . 

Remark: Note we have assumed that e > 0. The 
above system will fail to be ensemble controllable if 
e £ [— eo, £o]j as we cannot approximate an even func- 
tion /(e) = 9, with an odd degree polynomial. 

Remark The key idea in designing compensating 
pulse sequence is to synthesize higher order Lie brackets 
that raise the dispersion parameters to higher powers. 
The various powers of the dispersion parameter can be 
combined for compensation as explained above. The con- 
struction presented here is not the most efficient way of 
achieving a desired level of compensation. The construc- 
tion given here however presents in a transparent way 
the role of higher order lie bracketing. We now consider 
an example when there are more than one parameter in 
the system dynamics. 

Example 2 Now consider the system 

X = (e 1 u(t)fl x + e 2 v{t)n y )X. 

where e\ £ [1 — <5i, 1 + 8{\ and £2 £ [1 — #2> 1 + £2]) fc> r 
< 6 1 < 1 and < 82 < 1. The system is ensemble 
controllable with respect to dispersions e± and e 2 - 

The reasoning proceeds along the same lines as before 
except now we have two dispersions parameters that are 
independent. Let adx(Y) represent the lie bracket [X, Y] 
(similarly ad 2 x {Y) = [X, [X,Y]]). Consider the identity, 

ad 2 X(e 2 Q y )^(-lfe l 2k + 1 e 2 Q z . (2) 

for k = 0, 1, 2, . . . , n. 

We can now choose coefficients Ck such that 
Sfc c fe e i' c+1 approximates a constant function over the 
range of t\. As a result, we can generate the bracket di- 
rection 62^2. Now using the bracket directions 62^ and 
62^l y , and the construction in Example 1, we can further 
compensate the dispersion of €2 and steer the whole en- 
semble together to a desired point. Infact the final point 
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can be made to depend explicitly on e\ and e 2 by synthe- 
sizing the bracket directions J2ki( c ki e i k+1 e 2 l+1 )^z , an d 
J2ki(dki£i k £2 l+1 )^y The coefficients Cjy and d a can 
be now so chosen that we can approximate rotations 
exp(0(e 1; e 2 )Q x ) and exp(#(e 1; £2)^1,). Therefore we have 
ensemble controllability. 

Example 3: Phase dispersions cannot be com- 
pensated Consider an ensemble of Bloch equations 

X e = A{t)(cos(</)(t) + 9)Q X + sin(0(<) + 9)Q y )Xe, (3) 

where there is dispersion in the phase of the rf field. The 
system is not ensemble controllable with respect to the 
dispersion 9 S [#1,6*2]. 

Proof: The simplest way to see this is to make the 
change of co-ordinates Y$ — exp(—£l z 9)Xg. The resulting 
system then takes the form 

Y = A{t){cos{<t>{t))n x +sm{<f>(t))n y )Y e . 

Since all Yg see the same field, they have identical trajec- 
tories. As a result Xg cannot be simultaneously steered 
from (0, 0, 1) to (1, 0, 0). Lack of ensemble controllability 
can also be understood by looking at Lie brackets of the 
generators. Equation can be written as 

X e = {A{t) cos(0(t)).Bi + A{t) sm(4>(t))}X e , 

where the B\ = cos(9)fl x + s'm(9)fl y and B 2 — 
- sm(9)n x +cos(9)n y . Observe that S 3 = [B 1 ,B 2 ] = fi z . 
Therefore, all iterated brackets of B[s are linear in cos(6>) 
and sin^) and we cannot raise the dispersion parameters 
008(6*) and sinf^) to higher powers and therefore cannot 
compensate for the dispersion in 9. 

Example 4: Larmor dispersion in the presence 
of strong rf-field Now consider the Bloch equations 

X e = (cofl z + u(t)Cl x + v(t)Q, y )Xg. 

with dispersion in the Larmor frequencies. The system 
is ensemble controllable with respect to the dispersion 
parameter to. 

Note because of the assumption of strong fields, we can 
reverse the evolution of the drift term 

e-xjp(nQ x ) exp(ujQ z dt) exp(— tt£I x ) — exp(—ujfl z dt). (4) 

Now as before a maneuver 

exp(— ujflzVdt) exp(-~fL x Vdt) exp(ujU z Vdt) exp(n x Vdt) 

produces the bracket direction [u;f2 z , £l x ] = loVL v to lead- 
ing order. Similarly [u;f2 2 [a;fi z , ft x j] = — (u>) 2 Q x . Hence, 
we can generate higher brackets with even and odd pow- 
ers of w. To see that the system is ensemble control- 
lable consider the Lie bracket relation ad^Q^inU.x — 
{-l) n uj 2n VL x and ad^n^+iMy = {-l) n+1 uj 2n + 1 n x , we 
can synthesize an evolution exp(^ fe CkU k Q, x ) and sim- 
ilarly the evolution exp(^ fc dk0J k Q y )- The coefficients 
Ck and dk can be chosen to approximate Euler angles 



(a(uj), 7(w)) an d we therefore as in Theorem 1, have 
ensemble controllability. 

Remark Note if we have only one quadrature of the 
control field i.e., 

X = (ujfl z + u(t)Q x )X, 

then we can only synthesize the generator £l y with odd 
powers of u>, with uj £ [—B,B], Therefore an evolution 
of the form exp(f (ui)£l y ) cannot be approximated if / is 
an even function. 

Example 5: Dispersion in Coupling Strengths 
Consider two coupled qubits with Ising type interactions 
with dispersion in coupling strengths J. The interaction 
Hamiltonian H c — Jai z o~ 2z , with J G Jq[1 — S, 1 + 5], 
5 > 0. Although not necessary, for simplicity of exposi- 
tion, we assume that we can produce local unitary trans- 
formation on the qubits much faster than the evolution 
of couplings. We now show that it is possible to compen- 
sate for dispersion in J and generate any quantum logic 
with high fidelity 

By local transformations we can synthesize the effec- 
tive Hamiltonian 

Jo\ y a 2z = exp(ia lx ^)(Ja lz a 2z )exp(-ia lx ^). 

Now using B\ = —i2a\ y a 2z and B 2 = —i2o\ z o 2z as gen- 
erators we get [JBi[JBi, JB 2 ] = — J 3 B 2 . Now using a 
construction similar to one in example 1, we can syn- 
thesize the evolution exp(^ fc Ck J 2k+1 cri z a 2z ), where the 
coefficients Ck are chosen such that J^k c kJ 2k+1 ~ Jo over 
the range of dispersion of J. Hence we have compensate 
for dispersion in J. We also have ensemble controllability 
with respect to the parameter J. Let 

A(J) = exp(-i(a( J)o\ x o 2x + b(J)cr ly or 2y + c{J)o lz o 2z )). 

We can write an arbitrary two qubit gate with the de- 
pendency on J as 

U 2 (J) <g> Ui{J) A(J) V 2 (J) g> Vx(J). 

where U\,V\ and U 2 ,V 2 are local unitaries on qubits 1 
and 2 respectively. We can synthesize them with a ex- 
plicit dependence on J as follows. Using the commuta- 
tion relations of the type [— i J2ai y a\ z , — i J2o\ z o\ z \ — 
—iJ 2 ai x , we can synthesize generators —i(J 2 )ai x , 
-i(J 2 ) k a ly , ~i{J 2 ) k a 2x , -i(J 2 ) k <r 2y (k = 0,1,2,... ) 
and use these to synthesize U\(J), Vi(J), U 2 {J), V 2 (J). 

Remark Using similar ideas as above, it is possible to 
compensate for more general coupling tensor. Consider 
the coupling tensor 

aa lx a 2x + (3<J ly a 2y + ja lz a 2z . 

with dispersion in a, (3, 7. Now observe for U = 
cxp(-ina x ), and A = cxp(-i(aai x a 2x + (3ai y a 2y + 
lo lz a 2z )), 

UAU ] A = cxp{-ij2a lz a 2z ). 

So we only need to take care of the dispersion in 7 and 
the construction is similar to the one before. 
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III. ENSEMBLE CONTROLLABILITY OF THE 
BLOCH EQUATIONS WITH BOUNDED 
CONTROLS 

We consider again the system Q but now with 
bounded controls, so that we cannot produce rotations of 
the type exp(— £l x n) in arbitrarily small time as in equa- 
tion (0J . Nonetheless the system is still ensemble control- 
lable as shown below. Our construction initiall y f ollows 
the well known algorithm of Shinnar-Roux mill. We 
then show how this construction can be extended to show 
ensemble controllability with respect to larmor dispersion 
and rf-inhomogeneity in Bloch equations. The solution 
to the Bloch equation Q is a rotation 

X(T) = RX(0), 

where R € SO(3). We work with SU(2) representation 
of these rotations. Recall a rotation by angle <fi around 
the unit vector (n Xl n y ,n z ) has a SU(2) representation 
of the form 



U 



a -(3* 
(3 a* 



(5) 



where a and (3 are the Cayley-Klein parameters satisfying 

a = cos in, sin — , o 

2 2 

(3 = -i(n x + in y ) sin-, 



aa* + 13(3* = 1. 
The Bloch equation then takes the form 



(7) 

(8) 



U = -2 



u u — iv 
u + iv — to 



U. 



The rotation U is simply represented by its first column 



(also termed spinor representation) ip 







Wc first 



consider piecewise-constant controls u(t) and v(t). The 
net rotation under these controls can be represented as 
successive rotations 



with the initial condition 



a 




' 1 " 


Jo 








The duration 



At, over which the controls u and v are constant can be 
chosen small enough such that, the net rotation can be 
decomposed into two sequential rotations since 



Under this assumption, we can write the rotation Uk as 
a rotation around z-axis by an angle coAt followed by a 
rotation about the applied control fields by an angle 
in SU (2) representation 



C k 



c k 



z 1 / 2 
z' 1 ' 2 



where 



c k 


9k 

= cos T , 


s k 


2 




= A k At, 


9 k 


, -1 v k 

= tan — , 

u k 








= J u l + V k > 







(10) 



^, (11) 



Plugging (|10|l into J5J|, we get the recursion relation of 
the spinor 



ot k 



C k 

s k c, 



S k z- 
k z 



Olk-l 

0k-l 



Defining Pk = z k / 2 otk and Qk = z k l 2 f3k 1 the recursion 
may then be reduced to 



(12) 



" Pk ' 




'c k 


SI*' 1 ' 




' Pk-i 






_ s k 


Ckz- 1 




_ Qk-l 



with the initial condition 



" Po ' 




' 1 " 


Qo_ 








(13) 



Having the recursion (|12|) and the initial condition JT3J|, 
the spinor at the nth time step can be represented as the 
{n — l)-order polynomials in z (the parameter z encodes 
the dispersion parameter to. 



where Ua 



U = U n U n - 1 ...U 1 U , 

and dj, bj are the Cayley-Klein 



a, -b* 
b 3 a* 



parameters for the jth interval. Defining the multiplica- 
tion of the matrices Uj up to k by 



OL k -PI 




a k 


-n~ 




a 


— h* ' 
°0 


Jk a% _ 




b k 


a t _ 




bo 


a* 



the effect of the controls can then be calculated by prop- 
agating the spinor 



(9) 



a k 




a k 






a k -i 


_(3 k _ 




b k 


a t 







Note that 



Pn{z) 
Qn{z) 



E 

k=0 
n—1 

E 

k=0 



PkZ 



„-k 



qkz 



„-k 



\Pn(z)\ 2 + \Q n (z)f 



(14) 
(15) 

(16) 



which follows from (JSJ. The representation of rotation 
produced by the controls u(t) and v(t) has now been 
reduced from a product of n matrices in SU{2) to two 
(n — l)-order polynomials. The desired final states of an 



5 



ensemble of systems in Q, described by Cayley-Klein 
parameters, are two functions of z, and hence of uj. Wc 
can now design two polynomials P n (z) and Q n (z) such 
that we can approximate any desired smooth functions 
F a (z) and Fp(z) satisfying \F a{z) \ 2 + \Fp{z)\ 2 = 1, which 
characterizes the desired spinor we want as function of z. 
Now we can work backwards and compute the u' k s and 
v' k s that will produce P n (z) and Q n {z). Note by multi- 
plying both sides of 112(1 by the inverse of the rotation 
matrix we get 



(17) 



and the constraint of Hit) [I is still preserved. We have a 
backward recursion where we use the knowledge of co- 
efficients of Pk{z) and Qk{z) to compute Pfc_i (z) and 
Qk-i(z). This is the well known Shinnar Roux {Ta, lb"H 
algorithm. Because Pk-i(z) and Qk-i(z) are lower order 
polynomials, the leading term in Pk-i and the low-order 
term in Qk-i{z) must drop out 



" Pk-l ' 




C k P k - 


- S^Qk 


Qk-i 




_ (S k Pk - 


1- C k Q k )z 



CkPk,k-i + S k Qk,k-i — 0, 
—SkPk,o + CkQkfi = 0, 



(18) 
(19) 



where Pk, m denotes the coefficient of z~ m term in Pk(z). 
Observe that these two equations are equivalent as may 
be seen by expanding (|16fl as a polynomial, 



n— 1 m 



P W -EE [PiPm-i + QiQ*m-i] 



m=0 i=0 



and noting that all but the constant term are zero. The 
coefficient of z~( fe -i) in P{z) gives 



Pk,k-lPk0 + Qk,k-lQ 



k,0 



0. 



With this relation either equation (|18fl or (|19|) may be 
derived from the other. Choosing (|19|) and combining it 
with l|12[l . we get 



ik,0 



-le 



cos 



Pkfl 

This gives the rotation angle 
(hi- = 2 tan -1 



!k,0 



k,0 



(20) 



(21) 



Combining 12U|) and (J2TJ, we obtain the phase of the 
controls 



6i, — Z 



V P fc ,o / 

The controls Uk and Vk are then 



<t>k . a 



(22) 



4>k 
At 



cos 6» fc , 



These expressions for controls coupled with the inverse 
recursion in (|17[1 construct the piecewise constant con- 
trols Uk,Vk that generate polynomial approximations 
P n (z) and Qn(z) of the target function F a (z) and Fp(z). 

In particular, if we choose Q n (z) = — zsin^ and 
P n (z) — cos ^, we obtain a broadband rotation around x 
axis by angle <fi and similarly by choosing Q n (z) = sin^ 
and P n {z) = cos^, we obtain an approximation to a 
broadband rotation around y axis by angle <fi. 

If the amplitude of the controls is bounded, we can 
choose 4> small enough so that it can be achieved by small 
flip angles 4>k in equation (|22|l . Now we can concatenate 
these rotations to achieve a rotation with a bigger angle 
and thereby maintain the bounds on the control. 

Now we consider the case when there is also rf- 
inhomogeneity. If we produce a small flip angles <f> com- 
pensating for dispersion in u>, then dispersion e in the 
strength of the control u and v, results in Q n {z,e) ~ 
—ie l9k ^e and the ensemble executes an effective rotation 
exp(— (cos^fc)^ -fsin^fc)^)). Now using methods of 
example one, we can concatenate many such rotations to 
compensate for e. We now show ensemble controllability 
with respect to both dispersion in the natural frequency 
u> and strength of the rf-field. 

We again write the final rotation U £ SU (2) as 



U =U n U n - 1 ...U 1 U , 



where 

U k (z,e) = 



C k (e) 
S fc (e) 



C fc (e) 



zV2 o 
z-V2 



Note that the flip angle has a dependence on the pa- 
rameter e. We can now choose a desired Q n {z,e) = 
Y2=l <lk{t)z~ k and P n (z,e) = J^kZo Pk{()z~ k and find 
the flip angles (^fe(e), ^fe(e)) that creates these polynomi- 
als. Now we can use the results of Example 1 to find 
pulse sequences that will synthesize (0(e), </>(e)). This 
then establishes the ensemble controllability with respect 
to both uj and e. Such constructions can also be used to 
generate pattern pulses that selectively excite the Bloch 
equations with parameters lying in a given subset of lu — e 
space 8j. 

□ 

We now investigate the subject of ensemble control 
from a general control theory perspective. 



IV. ENSEMBLE CONTROLLABILITY 



Consider a family of control systems 
f s (x s ,u,t), 



dx s 
IT 



(23) 



indexed by the parameter vector s taking values in some 
compact set f2 C R d . The same control u(t) E R m is 
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being used to simultaneously steer this family of control 
systems. For such systems, we define the notion of 
ensemble controllability as following. 



Definition 1 The family of systems in (I23|) is called 
ensemble controllable, if there exists a control law u(t) 
such that starting from any initial state x s (0), the 
system can be steered to within a ball of radius e around 
the target state g s , i.e. || x s (T) — g s ||< e. Here x s (0) 
and x s (T) are interpreted as functions of the variable 
s at time and T and || • || denotes a desired norm on 
function space. In this paper we take |j • || to denote L2 
norm. The final time T may depend on e. 



Definition 2 An ensemble of systems is called point 
ensemble controllable if both the initial state x s (0) 
and the target state g s are constant functions and the 
ensemble of systems can be steered between the initial 
and final states as defined above. 

The key problem of interest in ensemble controllability 
is to characterize the properties of the ensemble such that 
if each system of the ensemble is controllable, the whole 
family is ensemble controllable. We begin with consid- 
ering linear systems x = Ax + Bu, where B is a n x m 
matrix and u € R m 

Theorem 1(A negative result:) An ensemble of 
linear control systems is not ensemble controllable if there 
is a variation in the control matrix B. 

Proof: Consider the family of systems 

dx s , . 

= {Ax s + B s u). 

By variation of constant formula 

x s (T) = exp(AT)x s (0)+exp(AT) [ exp(-AT)B s u(T)d,T. 

Jo 

Since x s (T) is linear in control, variation in B in general 
makes the system ensemble uncontrollable. 

Remark Observe B s u = Y^k=i u kb k s where b k are the 
columns of B. we can think of bk as constant vector 
fields that generate translations. Since bk all commute, 
their Lie brackets do not generate terms carrying higher 
powers of the dispersion parameters. 

Theorem Consider an ensemble (A, b) s , of linear sin- 
gle input controllable systems 



exists a similarity transformation T s such that 



x = A s x + ub s . 



(24) 



The system is ensemble controllable only if distinct A s 
have distinct eigenvalues and A s is full rank over the 
ensemble. 

Proof: Since each of the systems is controllable, there 



A. = t: 








-a s o 



1 





-a? 



(25) 



and 



b = Tl 



(26) 



where af are the coefficients of the characteristic poly- 
nomial of A s and show dispersion over the ensemble. 
First observation is that for the ensemble to be con- 
trollable, we must have T s to be a explicit function of 
(«Iq, af, Oq, • • • , a®_ x ). Suppose two distinct T\ and T2 
correspond to the same (ao, a±, 02, . . . , a„_i), then start- 
ing from x — 0, for equation (|25|l . the final points 
x\ and X2 for the two systems at any time t are re- 
lated by x\(t) — TiT2 1 X2(t). Hence there exists no 
control that steers the two systems to the same point 
starting from 0. Therefore T s should be a function of 
(ciq, af , a|, ■ • • , Q-n-i)) i- e - no two distinct A have the 
same characteristic polynomial. We now show the ne- 
cessity of A to be full rank. Observe if A is not full 
rank, i.e., ag = 0, then A k b for k > always lies in 
the subspace H — span{A l b}"^ . Therefore the en- 
semble starting from x = 0, cannot be driven to a final 
state /(af , a|, . . . , af l _ 1 )b, where / is an arbitrary smooth 
scalar function of of (af , a|, . . . , a^_i). 

Remark As mentioned before the key idea in design- 
ing compensating pulse sequence is to synthesize higher 
order Lie brackets that raise the dispersion parameters 
to higher powers. As a results nilpotent control systems 
are not ensemble controllable as we cannot generate a de- 
sired higher power of the dispersion parameter. Consider 
the control system 



= '^2u t (t)e l g l (x). 



If the Lie algebra generated by gi is nilpotent, the system 
is not ensemble controllable. 

Example Consider the well studied nonholonomic in- 
tegrator 



d 
dt 



The system is not ensemble controllable with respect to 
the parameter e. The control vector fields g\ = — 

and <?2 = — x l^~ generate a nilpotent algebra 



X\ 




1 




' " 


X2 


= me 





+ u 2 e 


1 


. X3 . 




_ ~x 2 _ 




_ Xi _ 



Q x 
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the Heisenberg algebra. Observe that [egi, ec/2] = 2e 2 g§^-, 
which commutes with everything else. 

Remark Let G be a semi-simple Lie group and let 
IeG. Consider the control system 

X = (J2uieiBi)X, 

i 

such that Ci have dispersion in the range [1 — 6i, l + £,],for 
Si > 0. If Lie algebra generated by {Bi} spans the tan- 
gent space of G, then the system is ensemble controllable. 
The construction is similar in spirit to example 2. 

In this paper, we have tried to motivate the study of 
problems involving control of ensembles of dynamical sys- 



tems with dispersion in the parameters. Such problems 
arise naturally in areas of coherent spectroscopy, quan- 
tum information processing and control of quantum sys- 
tems in general. We have tried to make explicit the role 
of noncommutativity as a key aspect of the dynamics 
that makes design of a compensating control signal pos- 
sible. We note again that the constructions given in this 
paper donot provide the most efficient schemes for com- 
pensation, yet they illustrate the main ideas of the paper 
in a transparent way. The subject of ensemble control 
appears to be very rich and we anticipate that many in- 
teresting results of both theoretical and practical impor- 
tance will arise from a systematic study in this field. 
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